The contribution to the antimatter flux from individual dark matter substructures 
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The local antimatter fluxes induced by an individual dark matter (DM) substructure can be 
significantly dependent on the proper motion of the source. We derive analytic solutions to the 
propagation equation for time-dependent positron and antiproton primary sources, finding that the 
static limit is a fair approximation only for very high energy particles and nearby sources. We 
ON ■ discuss weakly interacting massive particle (WIMP) models fitting the PAMELA positron excess 

and the FERMI all-electron data. We show that, for a single non-static DM point-source, one 
cannot extract from the data, in a unique way, model independent particle physics observables, such 
cn : as the WIMP mass, the pair annihilation cross section, and the annihilation yield. The gamma-ray 

emission associated to WIMP models inducing a significant local flux of positrons or antiprotons is 
found to be compatible with EGRET measurements, but it can be definitely singled out with the 
FERMI-LAT telescope. 
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^ ! I. INTRODUCTION 

O ' Tlx- rwiil measurements of llie local cosmic-ray electron and positron Muxes have stimulated a lol of" inloresl in 

the existence of a nearby positron source in addition to the secondary positron component due to the interaction 
of primary species with the interstellar medium along propagation. The very recent high-statistic measurement by 
FERMI [2| of the all-electron flux (namely the flux of electrons plus positrons without charge discrimination) , while 
not confirming previous hints of an anomalous peak by ATIC Q and PPB-BETS [J], has found a spectral index 
sensibly harder than the one inferred from past extrapolations for the electron spectrum based on lower energy (and 
less accurate) data. Such a picture leaves room for a substantial contribution from unconventional (or previously 
disregarded) lepton sources, and it is very tempting to consider the possibility that a large component of the positron 
■ and electron flux is provided by the pair annihilation of dark matter (DM) particles, possibly forming the DM halo 
of the Milky Way (for early proposals in this direction see, e.g., [1, H, 0|; for reviews on indirect DM detection see, 
e.g., [1, [l(|). The latest fits of the data with such a component, see, e.g., G3 confirm that this picture is 
viable. Two general features of the DM-induced source were well-known even from early analyses: i) to match the 
level of the positron background, the local pair annihilation rate needs to be much larger than the one expected for 
thermal DM weakly interacting massive particles (WIMPs), within standard assumptions for the Universe thermal 
history and the local DM distribution; ii) having normalized the annihilation rate to the observed positron flux, 
there are very stringent bounds on the WIMP model from measurements of the local antiproton flux, disfavouring 
annihilation modes giving rise to a hadronic yield, and favouring the leptonic channels. DM models fulfilling these 
requirements have been proposed recently, see, e.g., [15, [T3|. For what concerns the first requirement, they mainly 
focus on a mechanism to account for a mismatch between the thermally averaged annihilation cross section at the 
freeze-out in the early Universe and a much larger annihilation cross section in the halo today. Alternatively, one 
could invoke an enhancement in the positron signal based on the presence of a local population of dense dark matter 
substructures, with the pair annihilation rate being large because the average of the number density of DM pairs 
is much greater than one half of the square of the mean DM number density (i.e., in terms of the local DM halo 
density p, (p 2 ) 3> (p) 2 )- When considering an average effect within many realizations of the Milky Way substructure 
population as extrapolated from current ACDM numerical N-body simulations, the mean enhancements in the local 
cosmic-ray fluxes are typically very modest, possibly below a factor of few [I3,[l(|. Large effects, at the level needed 
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to account for PAMELA and FERMI data, have been claimed instead in connection to a (few) single, very dense, 
nearby substructure(s) 0, El, El, H3| ; the price to pay in this case is that one has to refer to a configuration with a 
very small realization probability according to the N-body simulations, or to rely on a subhalo picture which is less 
standard. 

Both in discussing average effects from a full subhalo population and in tracing the effect of individual substructures, 
the approach of the recent papers and the vast majority of papers in the literature has been to ignore the fact that 
one is dealing with a system which is not static, and the emission and propagation of charged particles has been 
treated in the steady state limit. The distribution of substructures in the Galaxy is not rotationally supported; their 
typical velocity can be estimated from the total mass density profile of the Milky Way. Assuming for simplicity 
spherical symmetry for the Galaxy and an isothermal sphere for the subhalo phase-space distribution function, the 
velocity dispersion for such a distribution is simply equal to y/3/2 times the value of the circular velocity [2l|, i.e., 
assuming 250 km s^ 1 for the local rotational speed (22|, about 300 km s _1 . An object moving at such speed on an 
orbit perpendicular to the Galactic plane has a time of crossing of the diffusive halo region for cosmic-rays, say a 
cylinder with a 4 kpc half-height, of about 10 15 s. Such value is comparable to the typical confinement time for cosmic 
rays as estimated in the simplified "leaky-box" propagation models, and to the energy loss timescale for electrons 
and positrons [231 ]. This is clearly just a qualitative argument (actually not even referring to the most appropriate 
quantities, see the discussion below) to illustrate that the effect of proper motion of DM substructure can be relevant. 
Indeed, we will show that there is the possibility that local antimatter measurements reflect a transient due to DM 
annihilations in a subhalo, rather than a source to be modeled in the steady state limit. 

The paper discusses the main features of the local antimatter fluxes resulting from individual DM substructures, 
taking proper motion into account. Analytic solutions for the propagation equation, as appropriate for positrons and 
antiprotons/antideuterons, are presented taking into account the most relevant terms, referring to a simplified descrip- 
tion of the diffusion region, the interstellar medium and radiation field (these are the same kind of assumptions which 
are usually implemented for DM studies in the steady state limit, as well as, most often, to discuss electron/positron 
fluxes from astrophysical sources). We show results in a few benchmark configurations for the propagation parameters 
and for the orbit of the DM source. We introduce a few sample possible interpretations for the PAMELA and FERMI 
positron/electron data, illustrating how sensitive one becomes to different assumptions. In particular, we show that, 
contrary to the picture extensively discussed in recent analyses, it is no longer true that one can extract from the 
data, in a unique way, model independent particle physics obscrvables, such as the DM mass, the pair annihilation 
cross section and the annihilation channel. We also consider the gamma-ray signals associated to this scenario and 
compare with current limits as well as with the detection prospects in the upcoming future. 

The paper is organized as follows. In Section [Til we present the description of the particle propagation in the 
Galaxy. We discuss point-like DM substructures as sources of positrons and compare to PAMELA and FERMI data 
in Section ITlTl and Hv| respectively The contribution to the antiproton flux is presented in Section [V] In Section IVTl 
we compute other detectable extra features of the electron/positron flux, such as the associated radiative emission 
and the dipole anisotropy of the spectrum. Section IVlTl concludes. Details on the solution of the transport equation 
for positrons and antiprotons are reported in the Appendix. 



II. THE COSMIC-RAY PROPAGATION MODEL 



The random walk of charged cosmic-rays in the turbulent and regular component of the Galactic magnetic fields is 
usually treated through an effective description. Defining n the number density per unit total particle momentum of 
a given particle species (i.e., n{p)dp is the number of par ticles in the momentum interval (p,p + dp)), the propagation 
equation can be casted in the general form (see, e.g. [241] V 
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Here Q is the source term, including primary, spallation and decay contributions, D xx is the spatial diffusion coefficient, 
v c is the convection velocity, diffusive reacceleration is described as diffusion in momentum space and is determined 
by the coefficient D pp , p = dp/dt is the momentum gain or loss rate, tj is the time scale for loss by fragmentation, 
and T r is the time scale for radioactive decay. The problem is usually solved for stationary sources and assuming n has 
reached equilibrium (i.e. setting the left-hand side of the equation to zero), through a fully numerical integration of 
the general model, see, e.g., the Galprop 25] and the Dragon [26[ packages, or implementing (a chain of) semi-analytic 
solutions valid within a set of simplifying assumptions, see, e.g., [27l [2q|. In our analysis we will follow the second 
route, implementing however solutions of the diffusion equations valid for positron or antiproton primary sources 
which are non-stationary. As commonly done, we will assume a spatially constant diffusion coefficient and introduce 
a spatial average for the positron/electron energy loss rate; we will also neglect both convection and reacceleration, 
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sketching the relative energy dependent effects through an appropriate choice of the momentum scaling of the spatial 
diffusion coefficient. This very simplified scheme, which is flexible enough to introduce many study cases without 
neglecting any of the physical effects we wish to discuss, can actually be sufficient for a fair description of some of 
the key observables in cosmic-ray physics. E.g., Ref. [2!| introduces, for the Galprop numerical package, the case of 
a standard diffusive model with spatial diffusion coefficient of the form: 



£>(p)=/JA>l — j , (2) 

with P being the particle velocity in unit of the speed of light, R being its rigidity, and with the following parameter 
choice: D — 2.5 • 10 28 cm 2 s^ 1 , Rq = 4GV, 5 = 0.6 for R > R and S = for R < R (having neglected D pp , from 
now on we will simply label D xx as D). The diffusion region is treated as a cylinder extending from +hh to — hh 
in the vertical direction, with standard primary sources in a thin layer around z = 0, and up to Rh in the radial 
direction; in the example of Ref. 29], hh = 4 kpc and Rh = 30 kpc. Parameters are tuned to reproduce observational 
constraints, and in particular the relative abundance of secondary to primary components. Indeed, running the public 
available version of Galprop within this setup, we find a fairly good fit to the Boron over Carbon ratio data (reduced 
X 2 = 1.23 for R > 4 GV, considering the B/C measurements at high energy by ATIC [H, CREAM [p], HEA03 [H, 
and CRN [33], and having assumed a spectral index for primary nuclei of 2.35 and 2.1 for, respectively R < 40 GV 
and R > 40 GV), and to the antiproton over proton ratio as recently measured by PAMELA. We label this parameter 
choice as "model A", and take it as reference benchmark case. A comprehensive discussion of the dependence of 
our results on the propagation model is clearly beyond the scope of this paper; the main trends we will present are 
essentially insensitive to slight readjustments in the parameter space. One of the most relevant parameters for our 
discussion is the vertical size of the propagation region hh- It is well known that the local measurements of secondary 
to primary ratios are mostly sensitive to D^/hh rather than to each of the two parameters. The degeneracy would 
be broken by local measurements of the so-called "radioactive clocks", namely unstable secondaries, such as 10 Be 
as compared to 9 Be. Such data are however not very accurate at present. We consider two extreme setups: a thin 
halo model with ft,/, = 1 kpc and a thick model with hh — 10 kpc; Dq is rescaled to, respectively, 0.56 and 4.6 in 
10 28 cm 2 s _1 as we find simulating this model with Galprop and refitting the B/C ratio (reduced \ 2 = 1.10 and 1.22; 
the thin halo model is labeled as "model B", while the thick one as "model C"). Within these models the local 
primary proton and electron fluxes are reproduced as well (for the electron injection spectrum we take a spectral 
index of 2.30 at high energy, i.e. above E — 4 GeV). The spectra obtained with Galprop in our reference cases will 
be used to as background estimates in the next sections; the local electron and positron flux are computed with the 
spatially-dependent energy loss terms following from the standard templates for the interstellar radiation field and the 
magnetic field profile as implemented in the code. When computing primary components from WIMP annihilations 
we will introduce instead the simplification of a spatially-constant energy loss term, referring to an average value 
valid in the local neighborhood. In the energy loss configuration "H", which we assume as standard reference, the 
synchrotron and inverse Compton energy loss terms are driven, respectively, by an average magnetic field in the 
diffusion region which is about B = 6/iG, and a mean background starlight density U = 0.75 eV cm -3 (34|, on top 
of the cosmic microwave background component. For comparison, we will also consider a template in which both 
quantities are sharply reduced, assuming B = 1 /iG and U = 0.4eVcm -3 (we label this energy loss configuration by 
"L"). 



III. POSITRONS FROM A DARK MATTER POINT SOURCE 



Consider the limit of a point- like dark matter substructure, entering the diffusion region at the point at the time 
ti and moving along an orbit r p (t) (e.g., r p (t) = fi + v s (t — U) if one can approximate such motion as a straight line 
trajectory with constant velocity v s ), with dark matter made of WIMP dark matter particles of mass M x , annihilating 
in pairs with annihilation rate (crv) and positron yield per annihilation dN e + /dp. The positron dark matter source Q 
at the position r, momentum p and time t, takes the form: 

Q(f,t,p) = 6 3 [?-f p (t)]^r, (3) 

where T, the total dark matter annihilation rate in the source, contains all terms not depending on spatial coordinates, 
momentum and time: 

T = (av) [ dr s = (av) V s . (4) 

J 2 Af 2 2 M| V ' 
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FIG. 1: Left Panel: For three sample values of the positron momentum at injection {pa = 100 GeV, 500 GeV and 1 TeV) and 
as a function of the positron momentum p e as measured locally, the proper motion scale Ad (black dashed curves) is compared 
to distance d m at which the positron Green function is maximal (thick solid curves) and the range of values for which the Green 
function is 10% of the maximum (thin solid curves). The propagation model AH is assumed, see the text for more details. 
Right Panel: The same as for the left panel, but for po = 100 GeV and the six combinations of diffusion and energy loss models. 



Here p s (r s ) is the dark matter density profile within the substructure, and the expression after the equivalence sign 
defines the annihilation volume V s , having normalized p s to the reference value po — 0.3 GeV cm~ 3 . The positron 
number density per unit momentum is given by: 

p rt fP™ dN e + 

n(r,p,t) = J dt a J dp G (r,t,p;r p (t ),t ,po) (5) 

The Green function G is given in Appendix [XJ neglecting boundary conditions it is approximately in the form: 



G ^ 7 r3/2[A(p, J30 )]3 eXp ' 



d(t ) 
Hp,po) 



2 ' 



6[(t-t )-AT(p,p )} (6) 



where we introduced the distance d = \f — r p (to)\ between the source and the observer at the time t$, the energy 
loss timescale At = J dp/\p(p)\, and diffusion length A, defined through A 2 = 4 Jj° dp D(p)/\p(p)\. Looking at this 
expression, we can guess that including the time dependence in the propagation equation is relevant whenever the 
variation of d within the time At, say Ad ~ v s ■ At, is larger or comparable to A. Consider the high energy limit, in 
which the energy loss term scales like p(p) oc p 2 and the diffusion coefficient like D(p) oc p s , with 5 ~ 0.3 — 0.7; the 
square of the ratio between Ad and A goes like: 

in - 2 ^3 2 oo 1-g 100°- 6 - 5 (f-^) 2 m 
A 2 Dm 1-0.6 p\+ s 1 ' 

where D\ and p\ are reference values for the diffusion coefficient and energy loss rate at 1 GeV, respectively, 
10 28 cm 2 s _1 and — 10~ 16 GeVs -1 , i>3oo the substructure velocity in units of 300 kms -1 , pioo is the positron mo- 
mentum in the equilibrium distribution in units of 100 GeV, and 1Z = p/po the ratio between such momentum and 
the momentum of the positron at the source. Being energy losses very effective at high energy, the time spent in the 
system before losing most of the energy by particles injected with very high momentum tends to be very small. The 
distance traveled by a DM subhalo in the same timescale is very small as well, and, in this case, proper motion can 
be safely neglected. An analogous picture occurs when considering cases with small difference between momentum at 
injection and momentum at equilibrium, and thus with short timescale associated to the electron/positron transport. 
One sees from Eq. [7] that Ad <§; A for 1Z — > or large p. On the other hand, we expect proper motion to be relevant 
at intermediate to low energies. This effect is illustrated also in Fig. [TJ in the left panel, for a few values of the 
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FIG. 2: Left Panel: Local positron flux spectral shape for a point-source composed by 500 GeV WIMPs annihilating into 
monochromatic e + e~, and moving along a vertical orbit intersecting the Galactic plane at a short distance from the observer. 
Each curve refer to a different time of observation, with the distance d of the source from the observer playing the role of time 
variable. This is useful since we can compare with the spectra one obtains for static sources placed at the same distances. 
The positron flux is multiplied by p 2 d? to match the approximate scaling in the static limit. Right Panel: For a give time of 
observation (namely, fixing d — 1 kpc) , predictions for the flux within our benchmark setups for the propagation model. 



momentum at the source po, we compute A as a function of p and, correspondingly, find the value of d mi namely, the 
distance from us along the direction perpendicular to the Galactic plane at which the exact Green function G rea ches 
its maximum (thick solid lines). Neglecting boundary conditions as in Eq. ([6]), we would simply find d m = \/3/2 A; in 
the plot the scaling of d m with p gets more rapid when the distance approaches the vertical boundary of the diffusion 
region at hh — 4 kpc (the propagation model "AH" is assumed in the plot). Also displayed are the range of distances 
within which the Green function is larger than one tenth of its maximum and our estimate for Ac?, as defined above. 
In the right panel, we plot the same quantities for a single injected momentum (p = 100 GeV), but looping over the 
three propagation models and the two energy loss configurations as specified in the previous Section. In both panels, 
comparing d m and the 10% range to Ad, one sees that there are large p intervals over which proper motion effects 
are expected to be large. 

To discuss the contribution to the local positron flux $ e from a single point source, we start considering the sample 
case of a source moving with constant velocity v s = 300 kms , on a trajectory which is perpendicular to the Galactic 
plane, intersecting the plane at a small distance from the Sun (say 10 pc, the precise value has little relevance for the 
discussion). To emphasizes propagation effects, we refer to a dark matter model, of given mass and total annihilation 
rate, with prompt annihilation into e + e~ pairs, i.e. with a monochromatic positron spectrum dN e + /dp ~ 6(p — M x ). 
In Fig. [5] we fix M x = 500 GeV and show the shape of the induced local positron fluxes in a few sample cases. In 
the left panel, each curve refers to a different time of observation, but rather than labeling them by time steps, we 
indicate the distance d of the source from the observer, along the trajectory and at the time when the positron flux is 
measured (dashed lines for an approaching source, solid line for a source which has passed nearby and now is moving 
away; for reference, in the example we are considering, a shift of 1 kpc in d corresponds to a time interval At ~ 10 s, 
i.e. about an order of magnitude lower than the typical "escape time" for particles of these energies as extrapolated 
in simple "leaky box" propagation models). For comparison, we also show the case of static sources (dotted lines), at 
a given distance from the Sun. From the Green function solution (again sufficiently away from the boundaries of the 
propagation region) we can guess that, in the static limit, the peak of the product p 2 ■ $ e scales like d~ 3 . This is indeed 
what we find in the plot, where the quantity d 3 p 2 <!> e has been normalized to 1 in arbitrary units for d — 10 pc; the 
figure shows also the departure from such scaling in case of a non-stationary source, with the result for approaching, 
static and departing sources essentially coinciding only in the limit of small distances and large energies. For moderate 
to large distances the height of the peaks increases (decreases) for a source which has passed by (is approaching). 
The spectral shapes are also very different, with the sharp cutoffs at low energy enforced by the mismatch between 
the time interval from positron emission to detection and the energy loss timescale, which increases at lower energies. 
Positron propagation is treated according to model "AH" as introduced above; such model has hh = 4 kpc and hence 
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FIG. 3: Local positron flux spectral shape for a point source made of 100 GeV (left panel) or 1 TeV (right panel) WIMPs 
annihilating into monochromatic e + e~ , in the benchmark propagation model AH. The two values of d refer to the distances 
which the source has from the observer at the time the flux is measured and after the source has passed in the vicinity of the 
observer; the set of curves displayed refer to a few sample orbits for the source (see the text for details). 



an approaching source at d = 3 kpc, along the vertical trajectory we are considering, has just entered the diffusion 
region and induces a local positron flux which is marginal compared to the flux from the same source after passing 
by and being on the other side of the trajectory at the same distance. Even after the source has left the diffusion 
region there is still a population of positrons which has been left behind contributing to the local flux, especially at 
low energies: here a larger d just reflects a longer time interval since the injection time and hence a more efficient 
degrading of the injected high energy positrons to small energies. In the right panel of Fig. [21 we take one of the 
cases plotted in the left panel, namely a source moving away from the observer and being at d = 1 kpc at the time 
of observation, and show instead the dependence of the spectral shape and the normalization of the measured flux 
on the benchmark models chosen for propagation and energy losses. Indeed the differences between the various cases 
are rather striking; this gives a feeling for the fact that there are large uncertainties related to propagation in our 
problem, which are not resolved tuning the model, as we did, to local measurements of secondary to primary nuclei 
ratios. Furthermore, there are potentially other relevant issues which we do not address in our simplified approach. 
E.g., we expect some dependence of the results on how boundary conditions are treated in the diffusion model: we 
are adopting here the standard approach of a spatially constant diffusion coefficient and free escape (i.e. the diffusion 
coefficient going to infinity) as a sharp transition at the vertical boundaries; this is the same approach followed, e.g., 
in Galprop. Considering a less sharp transition between the two regimes, with, e.g., an exponential vertical scaling 
of the diffusion coefficient as proposed by [2(| , we would obtain smoother cutoffs in the spectra of sources at large 
distances shown in the left panel of Fig. [2] A gradient in the energy loss term would also mildly affect the prediction 
for the local positron flux. 

Fig. [3] considers the cases for a lighter and heavier WIMP (respectively, 100 GeV and 1 TeV) and, referring again 
to the propagation model AH, discusses the dependence of the locally measured positron flux on the orbit of the 
source, having chosen two sample values of the present distance of the source from the observer (1 kpc and 3 kpc 
and source moving away from the observer). The thick solid lines correspond to the source moving along a vertical 
orbit intersecting the Galactic plane at 10 pc from the observer (i.e., the same orbit as in Fig. [2]), while thin solid 
lines are obtained in a few cases in which this scale is varied up to 1 kpc for d = 1 kpc, and to 2 kpc for d = 3 kpc, 
or inclining the orbit from vertical to a 45° incident angle; as it can be seen, such cases are essentially equivalent 
and show that what is actually relevant for the result is how much time the source spend within a distance from the 
observer corresponding to diffusion length A (with A depending on the emission and observation energies). Such time 
is about the same for all solid lines, while it changes if the source speed is changed. The dashed and dotted lines refer 
to the same orbit as for the thick solid line but with a source speed which is, respectively, twice as large and half of 
it, namely 600 km s _1 and 150 km s _1 versus 300 km s _1 . The dash-dotted lines refer instead to a circular orbit with 
velocity matching the observed circular velocity, i.e. about 250 km s^ 1 . For comparison spectral shapes in the static 
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FIG. 4: Two examples of fits of the PAMELA positron fraction (left panel) and of the sum of the electron and positron fluxes 
(right panel) with a component due to DM annihilations in a substructure. The monochromatic e + e~ final state of annihilation 
has been considered, as well as two sample values of the DM mass. The fit was performed assuming that the primary electron 
spectral index and normalization follows from the FERMI data. 



limit are also displayed. The picture emphasizes further the fact that the encounter with a dark matter point-source 
determines a transient in the local positron flux, and that the static approximation is roughly valid only for nearby 
sources, for which the energy of the measured positrons is not significantly different from the energy at emission, or 
for very energetic positrons. 

We have focussed the discussion on DM sources with a monochromatic positron spectrum, still it can be very 
simply extended to DM sources with a generic positron emission spectrum (which can be obviously thought as a 
superposition of line spectra for different masses and annihilation rates), along the same patterns regarding the 
dependence on distance and energy. 



IV. A NON-STATIC POINT SOURCE VERSUS PAMELA AND FERMI ELECTRON/POSITRON DATA 

As an application of the general discussion outlined in the previous Section, and to introduce a focus on a specific 
case despite the many ingredients and parameters involved in the problem, we consider the possibility of a major 
contribution from a single non-static DM point-source to the local electron/positron spectrum as recently measured 
by PAMELA Q and FERMI 0. 

For what concerns the background from ordinary astrophysical electrons and positrons, we will consider two possi- 
bilities. In the first scenario, the bulk of the "all-electron" spectrum measured by FERMI is due to primary electrons 
emitted in supernova cosmic-ray sources (with only a mild contamination, at the 10-20% level, from secondary elec- 
trons and positrons), and, hence, we can use this data-set to derive the electron spectral index at the sources. The 
second possibility is that the FERMI measurement has actually found an extra electron and, possibly, positron source, 
comparable to standard contributions at about 100 GeV or so, and dominating at high energy up to the cutoff found by 
HESS at about 1 TeV [35| : in this second scenario, we assume that the primary electron spectrum is softer and can be 
inferred from preliminary (less accurate) lower energy data on the electron-only spectrum presented by PAMELA 36] , 
while the extra contribution is dominated by a single DM source. 

Starting with the first background choice, we show in Fig. [4] (left panel) two examples of fits of the PAMELA 
positron fraction at high energy with a component due to DM annihilations in a subhalo. The fits have been obtained 
with sample values of the WIMP mass and a given annihilation channel, respectively 500 GeV and 1 TeV and the 
monochromatic e + e~ final state. We have also assumed that the source moves with constant velocity v s = 300kms _1 
on the trajectory perpendicular to the Galactic plane already introduced for Fig. [2j Looping over the propagation 
models discussed in the previous Sections, we have extracted the best-fit values for the source distance and the total 
annihilation rate in the source T. For a given propagation model, the background is inferred normalizing the proton 
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FIG. 5: Two further examples of fits of the PAMELA positron fraction (left panel) and of the sum of the electron and positron 
fluxes (right panel) with a component due to DM annihilations in subhalo. With respect to Fig. [4] we are now considering 
a different annihilation channel, i.e. t + t~, and that the primary electron spectral index is slightly harder than the one you 
would infer from the FERMI data. 



flux to the locally observed spectrum, making a prediction for the secondary electron and positron components using 
the Galprop package, and deriving the primary electron spectral index and normalization from the FERMI data. 
Including then the DM component, the fit to the PAMELA, FERMI and HESS data (the HESS datasets are rescaled 
within systematic uncertainties to match the normalization from FERMI) is performed allowing for a 20% variation 
in the normalization of the secondary spectra (reflecting various uncertainties on their modelling) and on the primary 
electrons, as well as a slight tilt in the primary electron spectral index; the fit disregards PAMELA datapoints below 
10 GeV, since the force-field method we are implementing to take into account solar modulation, with modulation 
parameter of 0.55 GV, is probably not sufficiently accurate, and a modulation which takes into account the charge 
sign should be implemented instead [l|. For a given WIMP mass, and for all propagation parameters except for the 
thin halo model (model B), this procedure produces configurations with moderate to small x 2 . They are restricted 
in a rather narrow range of allowed distances; however this range shifts with WIMP mass, choice of the propagation 
model and the velocity of the source along the orbit. The parameters for the two sample models obtained from the 
fit and plotted in Fig. [4] are given in the first two rows of Table |U 

Analogously, Fig.[5]shows two sample fits of the data in the case of a primary electron spectral index slightly softer, 
according to the picture mentioned at the beginning of this Section. Since we need to account for a positron/electron 
exotic contribution up to 1-2 TeV, DM WIMPs needs to be rather heavy (e.g., we take 5 TeV and 3 TeV), while 
the choice of the annihilation channel is again not critical (we chose to refer to the case of pair annihilation into the 
r + t~ final state). The results from the fit are reported in the last two rows of Table fl] A few issues should be 
stressed looking at this table. First, for all the four selected models, the \ 2 indicate fairly good fits, although not 
exceptionally good (one should also take into account that, especially for FERMI, errors in the data-set are correlated, 





M x 
GeV 


annihilation 
channel 


r 

10 36 s" 1 


v s 

kpc 3 


prop, 
model 


d 
kpc 


$ 7 (£ > 0.1 GeV) 
cm -2 s _1 


x 2 

(d.f-50) 


1 


1000 


e+/e- 


20.9 


5.3 • 10 5 


AH 


4.25 


1.2 • 10~ 8 


46.7 


2 


500 


e+/e~ 


73.4 


4.6 ■ 10 5 


CH 


6.25 


1.6 ■ 10" 8 


42.3 


3 


5000 


t+/t- 


1.9 


12.6 ■ 10 5 


AH 


1.54 


1.1 • 10" 8 


44.4 


4 


3000 


t+/t- 


2.4 


5.5 ■ 10 5 


CH 


1.43 


1.4 • 10~ 8 


60.9 



TABLE I: Sample DM models 



9 




, -15 -10 -5 5 10 15 20 

d [ kpc ] d [ kpc ] 



FIG. 6: Upper limits on the total annihilation rate F derived from the PAMELA data on the positron ratio and FERMI data 
on the all electron flux. A sample WIMP model is considered (M x = 500 GeV and monochromatic e + e~ final state), while we 
loop over propagation model configurations (left panel) and a few possibilities for the point source orbit (right panel). Results 
are compared to the level of P required to match either the brightest unidentified EGRET source, the faintest source detected 
by EGRET and the level of the 7-ray sensitivity for FERMI. 



while we have simply added in quadrature statistical and systematic errors); our interest is, however, to show the 
feasibility of the framework, rather than playing with all possible uncertainties in the background and the signal to 
improve the fits further (and indeed a slight readjustment of the background could lead to better fits in all cases). A 
second issue is that the values which we find for the total annihilation rate T are very large for all models, as one sees 
converting them to the annihilation volume V s under the hypothesis that the annihilation cross section is of the order 
av ~ 3T0~ 26 cm 3 s _1 (this is, roughly speaking, the level needed for a thermal relic WIMP to match the DM density in 
the Universe in case of standard assumptions for the Universe thermal history, and without invoking an enhancement 
in the cross section going from the freeze-out to the zero temperature limit, see, e.g., Q). Annihilation volumes of 
the order of 10 5 — 10 6 kpc 3 are much larger than typical values predicted for DM substructures in N-body simulations 
of hierarchical clustering for Milky- Way size DM halos, with a realization probability for a configuration containing 
these sources, supposing one can extrapolate from the results shown in Ref. [lj| for static sources, below few in 10~ 4 . 
On the other hand, in a different scenarios, such as the one in which the adiabatic formation of an intermediate mass 
black hole drives a sharp enhancement of the dark matter density inside the hosting substructure [37| , the probability 
density of V s for these sources has a peak at about 10 6 kpc 3 and a tail extending to much larger values [15||. There is 
also the possibility that the reference value av ~ 3T0~ 26 cm 3 s -1 is a significant underestimate of the pair annihilation 
cross section, and in this case V s would be appropriately downscaled. 

Having shown that it is not manifestly implausible to interpret current positron/electron data within the single 
DM substructure scenario, the examples chosen for Figs. 0] and [5] question the standard methods applied to extract 
model-independent informations on the DM candidate: It is usually assumed that interpreting an excess as a DM 
signal gives an estimate of the DM particle mass when the energy threshold for the excess is detected, that the spectral 
shape of the excess determines the dominant annihilation channel, while the normalization is mainly an indicator of 
the level of the pair annihilation rate. Here, instead, the threshold is just setting a lower limit to the DM particle 
mass, since the energy at which the exotic component dies out depends also on the distance of the DM point-source. 
Regarding the shape of the spectrum, we have shown that this is a transient, keeping little memory of the spectrum 
at the source, and indeed, in the plot, the spectral shapes of the single source contribution and that from the smooth 
DM halo component are sensibly different in all examples. Finally, the level of the induced flux depends mainly on 
the product av ■ V s , and it is difficult to give any solid estimate for the annihilation volume V s . 

In Table Q] we quote the prediction for the gamma-ray flux, integrated above 100 MeV, induced by the DM 
point-source in the four selected models: in the first two cases, with only monochromatic e + /e~ as tree-level final 
state of annihilation, the flux is due to photons emitted as a final state radiation (FSR), with the rate and photon 
energy distribution for the \X ~ * e + e~7 which are estimated in terms of the lowest order process XX ~^ e + e~ and 
in the approximation for m e -C M x (see, e.g., [38| and references therein; an eventual model-dependent "internal 



10 




-15 -10 -5 5 10 15 20 



d [ kpc ] M x [ GeV ] 

FIG. 7: Upper limits on the total annihilation rate F derived from the PAMELA data on the positron ratio and FERMI data 
on the all electron flux. The left panel is specular to the right panel of Fig. [7J but for the r + r~ final state. In the right panel 
the mass scale for the DM candidate is varied and a few sample values of the distance for a source moving away from the 
observed along a vertical orbit are considered. Results are compared to the level of the F required to match either the brightest 
unidentified EGRET source, the faintest source detected by EGRET and the level of the 7-ray sensitivity for FERMI. 



bremsstrahlung" contribution [3^| is not included here). For the last two cases, the annihilation into t + /t~ gives rise 
to decay chains containing neutral pions, which in turn decay to two photons; this process is accounted for by linking 
to the Pythia Montecarlo simulations as provided by the DarkSUSY package j40|. The level of the integrated gamma- 
ray fluxes, slightly above 1 • 10~ 8 cm" 2 s -1 , is much lower than the one for the brightest unidentified EGRET source, 



i.e. about 7 • 10" 7 cm" 2 s" 1 , or the level of the faintest source detected by EGRET, i.e. about 6 • 10~ 8 cm -2 s _1 41 1 
(in a more detailed comparison one should consider also the fact that the energy spectra for the sources proposed 
here are sensibly softer than typical spectra for EGRET sources) , but well within the sensitivity for the FERMI LAT 
instrument, about 4 ■ 10~ 9 cm" 2 s _1 considering intermediate to high latitudes and a 2 year data taking period [42j 
(again we are not taking into account spectral features, and also the fact that the sensitivity we quote for a DM signal 
was computed based on a gamma-ray background as extrapolated from EGRET GeV measurements of the diffuse 
gamma-ray emission which are not being confirmed by FERMI and will have to be sensibly lowered, see, e.g., fill). 
This is a trend we see in general. 

In Fig. [SI we compute the upper limit to the total annihilation rate F from the PAMELA data on the positron 
ratio and from the FERMI data on the "all-electron" flux (limits are extracted separately for the two datasets, 
allowing for extra freedom in the background normalization and spectral index; the most stringent one is displayed) . 
These constraints are compared to the values of T corresponding to an integrated 7-ray luminosity at the level of the 
brightest unidentified EGRET source, of the faintest source detected by EGRET and of the FERMI 7-ray sensitivity. 
We have chosen a benchmark WIMP model, with M x = 500 GeV and annihilating into monochromatic e + e~ . We 
display the results for: i) the reference vertical orbit and different propagation models (left panel, the displayed 
cases are the same as in Fig. [2j see the relative discussion in the text; the values of d refer to the distance of the 
point source from the observer at the time when the positron flux is measured, after the source has passed in the 
neighborhood of the observer); ii) the reference propagation model AH and a few possibilities for the point-source 
orbit (right panel, the sample orbits are a subset of those considered in Fig. [31 negative d refers to an approaching 
source, positive d to a source moving away from the observer). In the right panel of Fig. [6l we also sketch the 
error that can be induced by estimating positron constraints assuming the stationary limit. Dotted curves refer to 
static point-sources and show (making the comparisons with the appropriate color coding) that the bounds, except 
for nearby sources, are systematically overestimated for approaching sources and greatly underestimated for sources 
moving away from the observer. Clearly, the gamma-ray luminosity of a point-source scales with 1/dr. The scaling 
with distance, or better with time, of the local positron flux depends on many different ingredients which are hard to 
extract from observations. The picture, however, is intuitively clear: unless the DM point-source is extremely close 
to the observer inducing a very large gamma-ray flux (at a level which can be already excluded by EGRET), it can 
generate a substantial contribution to the electron/positron flux measured by PAMELA and FERMI, without being 
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FIG. 8: Spectral shapes for the local antiproton flux due to a dark matter point source moving on a vertical orbit with constant 
velocity v s — 300 km s -1 and for different values of the distance d of the source at the time of observation. In the left panel, the 
source is getting closer to the observer. In the right panel, the source, after passing in the vicinity of the observer, is moving 
away. A WIMP mass of 500 GeV and pair annihilations into W + W~ have been assumed, as well as the propagation model A. 
For comparison, spectral shapes for a static source at the corresponding distances and for pair annihilations in the smooth DM 
halo are also shown. 



in conflict with present gamma-ray observations. In particular, we expect, for intermediate distances, the balance 
to lean towards the electron/positron side since there is a time interval during which electrons/positrons go through 
a diffusion transient while the gamma-ray intensity decreases as the inverse of time squared. When the diffusion 
transient is over and the local positron population becomes negligible, i.e. at large positive distances in the plot, the 
gamma-ray limit takes over again. On the other hand, if a substantial contribution to the positron flux measured by 
PAMELA is actually due to a DM point source, it is unlikely that such source will not be detected by FERMI in 
7-rays. This does not necessarily mean that it will be identified as a DM source. Indeed, FERMI will measure spectra 
only up to about 300 GeV and spectral features might not be easily identifiable, as well as the cross correlation with 
the positron flux is not anymore unambiguous, as we already stressed. 

These conclusions do not depend critically on the WIMP model, as we show in Fig. [7J In the left panel, we consider 
the same configurations of the right panel in Fig. [6l but now for a WIMP model annihilating into t + t~ . In this 
case, the 7-ray yield is enhanced (being the 7r°-decay channel open) and the electron/positron yield is softer; still, 
at intermediate distances, a picture inducing a sizable local flux of positrons is not in conflict with the EGRET 
observations, but again within the FERMI 7-ray sensitivity. In the right panel of Fig. [JJ we show the constraints 
on r obtained considering a few values for the source distance (still moving away from the observer on a vertical 
orbit with a velocity of 300 km s _1 ) and varying the WIMP mass. The interplay between electron/positron and 
gamma-ray observations is changing only slightly with mass, while it is again evident that the static approximation 
may be misleading. 



V. ANTIPROTONS FROM A DARK MATTER POINT SOURCE 



We focus now the discussion on the possibility that dark matter pair annihilations within a substructure give rise 
to an antiproton yield. This is possible if the final states of annihilation include weak gauge bosons or quarks, instead 
of leptons only, as we considered so far. Analogously to the positron case, in the limit of a point-like dark matter 
substructure, the antiproton source function takes the form: 

Q p (r,t,E) = S 3 [r-r p (t)]^T, (8) 



where dNp/dE is the differential antiproton yield per annihilation and T is the total dark matter annihilation rate 
in the source already introduced above. Solving the propagation equation with this source function, we find the 
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corresponding number density per unit energy: 



n(f,E,t) = ^-T I dt G(r,t,E;f p ,t ), 



(9) 



where the Green function G is given in Appendix [Bj As for positron, it is instructive to consider its approximate 
form by neglecting boundary conditions; assuming at first that also antiprotons annihilations in the gas disc along 
propagation play a minor role, G goes like: 
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7T^[X p (E,t,t )f 



■ exp 
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[X p (E,t,t )f 



(10) 



where the diffusion length depends explicitly on time: X P (E, t, to) = 1\J (t — to) D(E). Suppose for simplicity that the 
substructure is moving along a straight line trajectory with constant velocity v s , namely f p (to) — r t + v s (to — t); in 
this case G can be rewritten in the form: 
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(11) 



where we have introduced a "static" timescale t s = d 2 /4D, with dt 



r\ being the distance between the point- 



source and the observer at the time of observation (which would be, of course, the distance between the point-source 
and the observer at all times in the static limit), and a "dynamical" timescale td = AD/\v s \ 2 , and we have defined eg 
as the cosine of the angle between the vectors v s and f t — r (which is negative for a source approaching the observer 
and positive for those moving away) . We recognize that td and t s are the appropriate quantities to compare for a first 
guess on whether proper motion is relevant or not in estimating antiproton fluxes from point sources; in particular, 
Eq. (jlip can be integrated over time, giving: 



dt a G 



1 
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(12) 



In the limit for td S> t s , this expression correctly reduces to the Green function of the three dimensional Laplacian, 
while for td ~ 2t s we see that proper motion starts to become important. For typical values of the propagation 
parameters in the model, we find: 



2t„ 



8D 2 (E) 



<%v 2 



1.5 



2 rp2S IQ2S- 



D{ E 



10 



As 



(13) 



.10(1 



where E\q is antiproton energy in units of 10 GeV and dt,i is the source distance in units of 1 kpc. As for the 
positrons, we find the effect of proper motion to be negligible either in the limit of nearby sources or going to high 
energies, while we expect it to be significant in the other cases. 

Actually, Eq. (|10|) and (llip contain an oversimplification, since they were derived neglecting antiprotons annihila- 
tions during propagation, an effect which defines a third timescale that could also be relevant: depending on energy, 
the antiproton loss via annihilation can be much larger than the leakage from the boundaries of the diffusion region. 
When including this extra timescale, the expression for G becomes less transparent; we show instead results for a 
few sample cases obtained by implementing the exact Green function. In Fig. [SJ we plot the spectral shape for 
the local antiproton flux from a source moving along the reference trajectory, namely a path perpendicular to the 
Galactic plane and intersecting it at a short distance from the observer, with the source moving at constant velocity 
v s = 300 km s^ 1 . We have also chosen the propagation model A, and defined the dark matter model through a sample 
value for the WIMP mass, M x = 500 GeV, and assuming that the annihilation is dominantly into W + W~ pairs. As 
in the previous plots, d refers to the distance of the source from the observer at the time of observation and is used 
instead of a time variable to compare more easily to the static limit (dotted curves in the plot). There is a evident 
transient at small and intermediate energies, soon after the source has entered the diffusion region; the spectrum 
starts to fully match the static limit case only when the source arrives very close to the observer. Moving away from 
the observer, the scaling with distance is less severe than in the static limit and a contribution to the local antiproton 
flux persists much later than the time at which the source has left the propagation region (the vertical boundary is 
at hh = 4 kpc in model A). Since high energy antiprotons diffuse more efficiently, at late times the spectral shape 
becomes steeper and starts again to differ sensibly from the shape of the component due to pair annihilations in the 
smooth dark matter halo (dash-dotted line in the plot). Note that the latter, contrary to the positron case, traces 
rather closely the shape of static point-sources. 
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FIG. 9: The locally measured antiproton and positron fluxes as a function of the distance of the point source (negative distances 
label an approaching source, positive values a source moving away) normalized to the fluxes at a closest encounter distance 
dmin = 50 pc and to the scaling with distance of the companion induced gamma-ray flux, i.e. 1/d 2 . W + W~ annihilation 
channel is considered (giving rise at the same time to antiprotons, positrons and photons) and three sample values of the WIMP 
mass M x . The antiproton and positron fluxes are plotted at the energies E which are equal to M x /5 (left panel) and M x /50 
(right panel). 



The W + W~ annihilation channel is a copious source of positrons and gamma-rays (mainly from the production 
and decay of pions). In Fig. [5J we compare the scaling with distance (i.e. time) of the locally measured antiproton 
and positron fluxes at an energy corresponding to a fraction of the WIMP mass (1/5 and 1/50, respectively, in the left 
and right panels) for three sample values of M x . Both antimatter fluxes are normalized to the scaling with distance 
of the induced gamma-ray emission, which goes simply like 1/d 2 . The ratio between charged-particle fluxes to the 
gamma-ray flux increases rapidly around d = since the local number density of positrons and antiprotons decreases 
less severely when the distance from the source is smaller than the corresponding diffusion length. At negative d, i.e. 
for a source approaching the observer, one can see again the transients due to the fact that the source has just entered 
the diffusion region; these transients are essentially specular in the two cases. On the other hand, the contribution to 
the local antiproton population tends to survive much longer than for the positron component. When the source is 
moving away, the antiproton flux scales, at least at intermediate distances and except for very energetic antiprotons, 
like 1/d 2 ; positrons have instead much sharper transients, due to the efficient energy loss term. 

Fig.[5]has been obtained for the reference propagation model (AH, with the "H" labeling the energy loss configuration 
for positrons only) and the reference trajectory. As for the positrons, we expect the results for the antiproton flux to 
depend on the choice of the propagation parameters and the orbit for the point source. This is illustrated in Fig. 1101 
in the left panel, we show the scaling of the locally measured antiproton flux with distance considering two sample 
WIMP masses, two sample energies and the benchmark trajectory, while looping over the three reference propagation 
models. The main effect here is due to the increase in the diffusion coefficient going from model B to model A and then 
to model C, as well as to the boundary conditions, entering more critically for the thin-halo model B (hh = 1 kpc); 
note, however, that even for model B a sizable contribution to the local flux may be due to sources which are well 
outside the diffusion region and that would not be included in the static limit treatment. In the right panel of Fig. [TT)I 
we consider instead a few of the orbits that we have already introduced to discuss positron fluxes, with the pattern for 
source velocities which is analogous here; note that in the circular orbit case the transient following a close encounter 
can be sufficiently long lasting for the induced antiproton population to persist up to the next close encounter (what 
we plot is the "equilibrium" configuration after a couple of full orbits). 

Finally, in Fig. llli we summarize the picture by comparing, for a WIMP model with W + W~ annihilation channel, 
the constraints on the total annihilation rate which are set by the latest measurements by PAMELA of the antiproton 
to proton ratio [Sij], the electron/positron measurements by PAMELA and FERMI, as well as searches for point 
7-ray sources by EGRET. Projected limits for the FERMI 7-ray telescope are also shown. Analogously to the 
case of the smooth DM halo component, we find for a DM substructure contribution that if WIMP annihilates in 
a channel in which antiprotons are copiously emitted, the measured antiproton flux sets tighter constraints than 



14 



10 



■a 

l a. 

s 

I a. 



10 



10 



10 



10 



10 





! 

200 GcV, M / 50 


Model A 
Model H 


1 TeV, M / 50 

200 GcV, M / 5 

1 1 TeV, M / 5 




\ « 1 / d : 










' ' i '? f 
■' <* -'I /' 




'' • '7 /• 
, •' i .7 1 ' 
.Mi 

/ ; • f \ 1 




'- '• ■ III I i 





-10 



10 



15 20 

d [ kpc ] 



3 

I 2. 

4 
S 

I 2. 



10 



10 



10 



10 



10 



10 



Model A 




90° 


600 


200 GcV. M / 50 




90° 


300 




cue 


- 250 


200 GcV. M / 5 




90° 


150 


1 ^ 


« 1 /d : 








''v *N 













-20 -15 -10 



10 15 20 

d [ kpc ] 



FIG. 10: Scaling of the locally measured antiproton flux with the distance of the point source from the observer (negative 
distances label an approaching source, positive values a source moving away), for two WIMP masses M x (200 GeV and 1 TeV) 
and two sample energies E (equal to M x /5 and M x /50), and assuming pair annihilations into W + W~. In the left panel, 
the results refer to the standard vertical trajectory and to our three benchmark propagation models A, B & C (intermediate, 
thin-halo & thick- halo, see text for details). In the right panel, for the propagation model A, the color coding refers to different 
orbits, namely three vertical trajectories with velocities of the source changing from 300 to 600 or 150 km s _1 , and a circular 
orbit in the Galactic plane with velocity matching the local circular velocity 250 km s . 



electron/positron data, except possibly for extremely heavy DM candidates. This trend is reinforced going to larger 
distances. The information from EGRET is inconclusive, while FERMI is going to cover most, if not all, the parameter 
space currently probed by antiproton or positron searches. The displayed results refer to a source moving along the 
reference vertical trajectory, approaching or moving away from the observer; however, analogous conclusions hold for 
other configurations. The extension to other WIMP annihilation channels also leads to similar results. 



VI. EXTRA FEATURES OF POSITRON/ELECTRON EMISSION FOR A DM POINT SOURCE 
A. A gamma-ray component from inverse Compton emission on starlight 

Radiative emissions are unavoidably associated to electron/positron yields; in particular, the Inverse Compton (IC) 
radiation component due to the interaction of high energy electrons on optical starlight is peaked in the range between 
tens to hundreds of GeV and may be sizable for the intense dark-matter lepton emitters we have been considering. 

The IC photon emissivity is obtained by folding, at any give n point within the diffusion region, the IC emission 
power with the electron/positron number density, see, e.g., [451 ] : 

j IC {v,?,t) = 2 dEP IC (r,E,v)n(r,E,t) , (14) 

J m c 

where the IC power is given by: 

f dft 

P w (r,E,v) = chu de-^{e,r)a{e,v,E)v (15) 

where e is the energy of the target photons, dn 7 /de is their differential energy spectrum and a is the Klein-Nishina 
cross section. The spatial dependence and spectrum of the number density of starlight photons can be computed from 
photometric maps of the Galaxy; we will adopt the model implemented in the public release of the Galprop numerical 
package (2f|. On the other hand, the analytic solution introduced in Section HTH for computing the electron/positron 
number density n was obtained assuming a spatially constant electron/positron energy loss term, and, hence, also 
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FIG. 11: Upper limits on the total annihilation rate F derived from the PAMELA data on the antiproton fraction (solid lines), 
PAMELA data on the positron ratio and FERMI data on the all electron flux (dash-dotted curves). The source is moving 
along the reference vertical trajectory and is approaching the observer (left panel) or moving away from the observer (right 
panel). Results are compared to the level of the F required to match either the brightest unidentified EGRET source, the 
faintest source detected by EGRET and the level of the 7-ray sensitivity for FERMI. 



a mean value for the background starlight density, rather than its value as a function of the spatial coordinates. 
This should not be regarded as a severe inconsistency, since the IC on starlight is just one of the effects contributing 
to the energy loss term, which main uncertainty is related to synchrotron emission on background magnetic fields. 
Having computed the IC emissivity field as in Eq. 1141 the relative 7-ray flux for a local observer is obtained simply 
by summing contributions along the line of sight. 

An example is given in Fig. 1121 where we refer again to the model introduced at the beginning of the discussion, 
namely a point source composed by 500 GeV WIMPs annihilating into monochromatic electrons and positrons and 
moving on a vertical orbit intersecting the Galactic plane at a short distance from the observer, see Figs. [Hand [5] The 
propagation scenario considered is the AH model (recall that in this model the vertical scale height of the diffusion 
region is 4 kpc) and the time variable is again replaced by the distance from the observer. Having assumed a velocity 
of the point source of 300 km s _1 , we consider three sample cases, i.e. d = 1, 3, 5 kpc with the source that has passed 
close to the observer and is now moving away. The induced IC flux (solid lines) is shown for three different directions 
of observation, all at the longitude I = 180° (i.e. opposite to the Galactic center), with latitude b = 90° (the direction 
towards the point source in our example), b = 75° or b = 60°. We plot also the background level corresponding to the 
emission from Galactic cosmic-rays (dotted lines) computed with Galprop, and the DM-induced 7-ray component due 
to FSR at the source (which is obviously present only in the direction of the DM substructure). The normalization 
of the total annihilation rate, or of the annihilation volume, has been chosen by saturating the current limits from 
PAMELA or FERMI, see the results in Fig. [5] This explains the relative strength of the FSR components. One can 
see that the IC flux in the direction of the source can be as large as the FSR contribution (for d — 1 kpc) or even 
larger (for d — 3 kpc). Most notably, the spatial size of the IC emission is significant, inducing a signal larger than the 
background even tens of degrees away from the source. Consider, however, that the case displayed here is the most 
favorable, since the source is located in the portion of the sky with the faintest Galactic background, and we have not 
included the uncertain extragalactic component. For d = 5 kpc, i.e. for the source outside the diffusion region, the IC 
term drops dramatically. This is due to the fact that there is no fresh source of high energy electrons and positrons, 
and the ones injected along the transit of the source through the diffusion region have rather rapidly degraded their 
energy. Therefore, there is no efficient high energy source to up-scatter starlight photons to 7-ray frequency. The 
trend is probably exaggerated by the fact that we are assuming a sharp boundary for the propagation region, while 
a smoother boundary condition would be certainly more physical and would probably lead to a smoother transition 
between the cases, d = 3 kpc to the d = 5 kpc. 

From the sample case discussed, we can infer that it is rather likely that if the local positron flux receives a sizable 
contribution from a single point-source, the induced IC flux is a relevant target for the FERMI 7-ray telescope. A 
more thorough discussion about this point is delayed to an upcoming dedicated analysis. 
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FIG. 12: 7-ray spectra for a DM point-source composed by 500 GeV WIMPs annihilating into monochromatic electrons and 
positrons, and moving along the reference vertical trajectory. The solid lines refer to the IC emissivity and are computed for 
three different directions of observation labeled by the latitude b (the longitude is I — 180°); dashed lines refer to the FSR 
component in the direction of the source, while dotted lines show the Galactic 7-ray background. The DM signals are computed 
for three sample distances of the source, assuming that, for each distance, the normalization of the total annihilation rate in 
the source is the maximum allowed by current PAMELA and FERMI electron/positron data, see the limits derived in Fig. [5] 



B. Dipole anisotropy in the electron/positron spectrum 

As pointed out by recent analyses (see e.g., [46|, 113, S, [4^] ) , a single nearby point-like source (e.g., a pulsar), being 
the dominant local positron source at high energy, induce an anisotropy which could be at a detectable level. In this 
Section, we reconsider this possibility in the context of DM substructure. We take into account proper motion effects, 
which are not relevant for bursting-like sources as pulsars. We assume the dominant contribution to the anisotropy to 
be given by the dipole term. In order to detect an anisotropy along a certain direction at a good confidence level, the 
required number of events has to be very large. Therefore, although the anisotropy in the positron spectrum would 
be higher than in the total electron plus positron spectrum, we refer to the latter, being the detection prospects for 
FERMI much more promising than for PAMELA. 

If the anisotropy is dominated by the dipole term, the intensity at a given point as a function of the direction of 
observation has only one maximum. The degree of dipole anisotropy can be defined by the quantity 8 = (I m ax ~ 
Imin) I \Imax +Imin), where I m ax and I m in are the maximum and minimum e + + e~ intensity with respect to direction. 
Expanding the intensity in spherical harmonics up to the dipole term, we have: 1(8, <fi) ~ 1+8 /cos 8+sin.d (I^ 1 e^ Ic ^ + 
I^e 1 ^), where 8 is the angle with respect to the direction of maximal anisotropy, / = (I max + Imin)/2 ~ = 
0) - 1(6 = tt))/2, and Jf 1 and ij are angle-independent coefficients. Note that contrary to the case of a stationary 
point-source (with homogeneous and isotropic propagation parameters) , the ^-dependent terms of the dipole are not 
null. 

Although the direction of maximal anisotropy (chosen, e.g., as the z-axis) is, in general, unknown, we can define 
the particle flux F along such direction by integrating over all the possible directions the projection over z of the 
intensity, i.e., F ~ Jdfll cos 8 ~ 4/3tt8I, and by estimating it as F ~ £>|Vn|. The latter is obtained in the 
diffusion approximation [50] . but it is valid at first order (i.e., for small anisotropies) also when including energy 
losses. Moreover, we assume the anisotropies in directions orthogonal to z to be subdominant (i.e., d z n ~ Vn). By 
equating the two expressions for F, one gets [51] : 

. 3D\Vntot\ 3D |Vn DM | , 1ft . 

= = t r , (16) 

c ntot c (n DM +ncR) 
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FIG. 13: Dipole anisotropy of the electron plus positron spectrum as a function of the energy for the four sample DM scenarios 
in Table U 



where udm denotes the contribution to the e + + e~ number density induced by a DM substructure, hcr is the 
contribution from cosmic-rays, and n tot — njjM + ncR- 

Neglecting boundary conditions, |VnDAf| is given by the gradient of Eq. [5] 



\Vn DM (r,p, t)\ 



E 



dt 



p(po) 



Pip) 



dN P + Xi 



dp A(p,p ) : 



G(f,p;f p {t ),po] 



1/2 



(17) 



with po obtained from At = t — to and x, L — x,y,z for i — 1,2, 3, respectively. 

We plot in Fig. [13J the degree of dipole anisotropy for the benchmark DM scenarios in Table Q] Comparing 
models 3 and 4, describing sources at analogous distance and with analogous local spectrum, one can note that, as 
intuitive, an approaching source would induce a higher anisotropy than a moving-away source (both on a straight- 
line trajectory). On the other hand, for sources at moderate distance the mismatch is quite small and the picture is 
similar to the stationary case. The degree of anisotropy is strongly and inversely correlated to the diffusion length (see 
Eg. [T7|. Although the benchmark DM models 1 and 2 induce a very similar local spectrum (see Fig. 2]) and describe 
substructures moving away from us on a vertical trajectory at the same velocity, in the second case the electrons and 
positrons undergo a greater diffusion before reaching us, and the anisotropy is suppressed with respect to the first 
case. Note that the maximum of the curves in Fig. [13] occur at larger energy than the peaks in Figs. [4] and [5) being 
the diffusion and the related wash-out of anisotropics less effective at high energy. The degree of dipole anisotropy for 
the benchmark models in Fig. [13] can be observable by FERMI in few years of data taking. On the other hand, the 
predictions for the anisotropy associated to a DM substructure suffer of many theoretical uncertainties, as in the case 
of the spectrum, and, depending on the trajectory, the degree of anisotropy can be suppressed or slightly enhanced 
with respect to an analogous (i.e., at analogous distance from the observer) stationary point-like case. 



VII. CONCLUSIONS 



We have discussed the contribution to the local antimatter fluxes due to individual WIMP DM substructures, 
accounting for the substructure proper motion. We have derived analytic solutions to the propagation equation 
as appropriate for a time-dependent positron and antiproton primary sources, identifying the relevant quantities to 
discriminate whether proper motion effects are relevant or not. We found that, for both positrons and antiprotons, 
the static limit is a fair approximation only in the case of high energy particles and nearby sources, while it fails in 
all other situations. The discussion has involved a few benchmark DM candidates and and a few sample orbits for 
the substructure. 

As an application of the general discussion, we focused first on WIMPs annihilating into leptons only, and we derived 
sample fits of the PAMELA positron excess and FERMI all-electron data. The fits have been obtained starting from 
some arbitrary values of the WIMP mass, demonstrating that, for a single non-static DM point-source, it is no longer 
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true that one can extract from the data, in a unique way, as extensively done in recent analyses, model independent 
particle physics observables, such as the DM mass, the annihilation yield and the pair annihilation cross section. 
Indeed, the threshold of the local electron/positron spectrum cannot be used as a sure indicator of the DM particle 
mass, being possibly set by the distance of the DM substructure. The relation between the shapes of the injection 
spectrum and of the propagated spectrum is not straightforward for a non-static DM substructure, which behaves 
as a transient source. Moreover, the annihilation rate depends on the product of the pair annihilation cross section 
times the unknown annihilation volume of the substructure, rather than on av only. 

The PAMELA positron excess can be explained in such a scenario, provided that a large annihilation rate in the 
substructure is considered. This requires either a large annihilation volume (much larger than typical values predicted 
for DM substructures in N-body simulations, but viable for scenarios accounting for the adiabatic formation of an 
intermediate mass black hole in the substructure) and/or a pair annihilation cross section significantly enhanced with 
respect to the reference thermal value. 

We have then used PAMELA and FERMI electron/positron data to derive, under different configurations, limits 
to the total annihilation rate in the DM source. These limits have been compared to the bounds extracted from the 
7-ray luminosities associated to the same DM models. The general trend is that the 7-ray signal induced by a DM 
point-source, giving a sizable contribution to the local positron flux, is below the level of the brightest unidentified 
EGRET sources, but well above the FERMI 7-ray telescope sensitivity 

We have also discussed WIMP models giving raise to an antiproton yield through pair annihilations in the sub- 
structure, as in the case with W + W~ as final state of annihilation. Analogously to the electron/positron case, single 
DM sources inducing a significant contribution to the local antiproton flux can be detected by FERMI in 7-rays (on 
the other hand, as for the smooth DM halo component, a positron/electron flux at the level of PAMELA and FERMI 
data can be hardly obtained in such WIMP models without violating antiproton limits). 

Finally, we sketched two further features in connection to positron/electron emission from a DM point-source, 
namely the spectral and angular shape of 7-ray flux induced by inverse Compton emission on starlight, and the dipole 
anisotropy in the electron/positron spectrum. We have shown that, potentially, both of them could be used to test 
the scenario proposed here. 
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APPENDIX A: SOLUTION TO THE PROPAGATION EQUATION FOR POSITRONS 

We present here a solution to Eq. [T] which is valid for a generic non-stationary primary source of electrons and 
positrons in a disc-like galaxy and with the following simplifying assumptions: We consider a model with diffusion 
coefficient and energy loss term that are a function of momentum only, while they are assumed to be spatially constant 
within the cylindrical propagation volume; we impose free escape at the vertical boundary of the diffusion region (i.e. 
that the equilibrium number density vanishes at z = ±/i^) and neglect the radial boundary (i.e. in the limit of radial 
size of the diffusion region Rh being much larger than hh). Reacceleration and convection are also neglected. In 
Cartesian coordinates Eq. |T]) becomes: 

dn{T Qp t] = Q(r,p, t) + D(p)An - |- (pn) (Al) 

The solution described here is a generalization of the one derived in Ref [52j to the case of diffusion coefficient D and 
momentum gain/loss term p which are generic function of the particle momentum p; the notation we use is hence 
analogous to that of [521 ]. We start by introducing N = —pn, and solve the equation for N rather than for n, using 
the Green function method and neglecting at first the vertical boundary condition: 

- ~ + -AG ~ = 5{r- r )S(t - t )S(p - Po ) . (A2) 
P at p op 



Going to the Fourier transforms in the variables r and t, we find: 

G = e £*>+**tos(p- Po ). (A3) 



-l{- lL o + D\k\ 2 ) - ^- 
p op 
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The solution of this equation is easily derived replacing the electron/positron momentum p with a double change of 
variables: 



dp 
Pip) 



and A 



f Po D(p) 
dfD(r) =-4 / dp- yl > 



Pip) 



(A4) 



and considering first the case p ^ po- The physically relevant solution vanishes for p > po since we will not include 
energy gain effects, but only energy losses (i.e., the term p(p) in the equations above is negative); we find: 

G = G e[^ AT - |£|2A2/4 ] ; (A5) 
where we defined At = r(p) — r(po) an d the constant Gq is obtained integrating Eq. (|A3[) around p = po: 



G 



(A6) 



Taking now the inverse Fourier transform, we find the Green function: 

exp 



Gfree(r,t,p;r ,t ,p ) 



A 2 



^3/2 A 3 



5 ((t- 1 ) -At) 



(A7) 



where it is now manifest to recognize A as a diffusion length scale and Ar as an energy loss timescale; the subscript 
"free" was added to the Green function since we have not included the vertical boundary condition yet. This is 
implemented with the image charge method, i.e. having defined z n = (— V) n z + 2nhh, and having factorized the Green 
function as G zb = G zb _ u ■ S ((t - t ) - Ar), G zb ^ ti is given by: 



G zb -u{r,p;f ,p ) 



exp 



(x-x ) +{y-yo)' 
F 



+ 00 



^3/2 A 3 



exp 



(Z n - Z ) 2 

A 2 



(A8) 



it is easy to verify that G zb (x, y, z = ±hh,t; ro, to) = 0. The density per unit energy is then: 



n(r,p, t) 



+h h 



-hh 



dx I dy / dz / dt / dp G zb (r,t,p;r ,t ,p ) Q(r ,t ,p ) . (A9) 



In case the source Q is point-like, such as the point-like dark matter substructure in Eq. ([3]), entering the diffusion 
region at the time t{, the density per unit energy becomes the expression in Eq. ([5]), or equivalently takes the form 
(which is the one we use for numerical calculations): 



n(f,p, t) —T dto 



p(Po) 



Pip) 



dN e , 
dpo 



■G zb -ti(r,p;r P (to),po) , 



(AI0) 



where po is obtained from the numerical inversion of the relation r{p) — t(pq) = t — to. Finally, we just remark 
that in case of a stationary diffuse electron/positron source, n(r, p) reduces to the familiar expression derived in the 
literature j53| (again generalized to D and p which are arbitrary functions of p) : 



n(f,p) 



Pip) 



+h h 

dx I dy I dz \ dp G z6 _ tl (f,p; r ,p ) Qir ,Po) , 

-h h 



(All) 



i.e. for a stationary axisymmetric source (i.e. for Q(ro,po) = QiRo, zo,po) with Rq the radial variable in a cylindrically 
symmetric system): 



n(r,p) 



dpo 
\Hp)\ Jo 



dRo Io 



2RR \ 2R ° ex P 



(R-RaV 
X 2 



+h, 

■hh 



+ OC 

dzo /J exp 



(z n - zof 



A 2 



Q(i?0, Z ,Po) , 

(A12) 



where we introduced the function Io{x) = cxp(— x) Ioix), with Io{x) the modified Bcsscl function of first kind. 
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APPENDIX B: SOLUTION TO THE PROPAGATION EQUATION FOR ANTIPROTONS 



Analogously to the positron case, we sketch here a solution to Eq.[T] which can be applied to a generic non-stationary 
primary source of antiprotons. We assume again that the propagation region is a cylinder and that the diffusion 
coefficient is spatially constant within it; we impose free escape as vertical boundary condition, while neglecting the 
boundaries in the radial direction. We still do not include reacceleration and convection, as well as energy losses give 
a negligible effect for antiprotons. We include instead the antiproton absorption effect, assuming it is mainly due to 
a thin gas disk, with gas target density constant in space and localized at z = 0. In Cartesian coordinates Eq. (fT]) 
becomes: 



dn(r, E, t) 

Ft 



Q(r, E, t) + D(E)An - P{E)5(z)n 



(Bl) 



with the term accounting for p losses due to collisions with the interstellar medium being P(E) — 2h g nHv(E)a p (E) , 
where h g is half of the gas disc height, and the gas disk is assumed to be made of hydrogen with constant density n H . 
The solution sketched here is analogous to the one presented in Ref. [54[ in a different context. It is derived through 
the Green function method, neglecting at first the vertical boundary condition: 



D(E) 



d 2 G d 2 G d 2 G 
+ 
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The following steps are to take the Fourier transforms in the x and y directions, i.e. G = dx J_°° dy Ge lk " x+lkvV 



D(E) 



d 2 G 

dz 2 



D(E)K 2 G — P(E)S(z)G 



dG 
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with K 2 = K 



k 2 , and to make a Laplace transformation in the t variable, i.e. G — J °° dtGe 
f) 2 C 

D—, - (DK 2 + s)G - P(E)S(z)G = - e ik * x °+ ik yV°- st °S(z - z ) 
dz z 

Away from z = and z = zq, the equation takes the form: 

d 2 G 



dz 2 



lo z G = 



(B3) 



(B4) 



(B5) 



with to = \J K 2 + s/D; the three solutions in the three distinct regions (i.e. the region between and zo, and those 
below and above it, setting the appropriate boundary conditions at z = ±oo) are matched to find G imposing its 
continuity and that (as one derives integrating Eq. (|B4j) around z = zq and z = 0): 
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One finds G = g- e ik * Xo+ih yy°- sto / D , where 

exp + \zq\)} \sign(z) - sign(z )\ 



fJ 
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2uj 2uj (P/D + 2uo) 

This expression can be inverted analytically to recover G. One can take first the inverse Laplace transform and find: 

q = e iKx +ik v y -{kl+kl)\l/i t . tQ ) ( B g) 

with: 
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where we defined: X p (E,t,t ) = 2yJ(t - t ) D(E) and L(E) = AD{E)/P(E). Finally, taking the inverse Fourier 
transform we find the Green function: 



Gf ree (f, t; fb, to) = —77 exp 
7rAr 



(x - x ) 2 + (y - y ) 2 



g(z,t;z ,to) 



(BIO) 



where the subscript "free" was added since we did not include the vertical boundary condition yet. As for the positron 
case, this is implemented with the image charge method, i.e. having defined z n = (—l) n z + 2nhh, the Green function 
is: 



G zb (r, t: f a , to) = -—x exp 

7TA p 



(x-x ) 2 + (y-yo) 2 



+00 



i" 1 )" 9(z n ,t;z ,t ) . 



(Bll) 



Note that the expression we found is analogous to Eq. lA8l for positrons; actually switching off antiproton annihilations 
(i.e. setting P(E) = 0) the two expressions are formally identical, except that the definition of diffusion length here 
is different from the diffusion length introduced for positrons, with time appearing explicitly in X p for antiprotons 
rather then implicitly through the momentum loss as in the A for positrons. The density per unit energy is then: 
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